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Abstract 



We develop series expansions for the ground state properties of the Hub- 
bard model, by introducing an Ising anisotropy into the Hamiltonian. For the 
two-dimensional (2D) square lattice half-filled Hubbard model, the ground 
state energy, local moment, sublattice magnetization, uniform magnetic sus- 
ceptibility and spin stiffness are calculated as a function of U /t, where U is the 
Coulomb constant and t is the hopping parameter. Magnetic susceptibility 
data indicate a crossover around 17 ~ 4 between spin density wave antifer- 
romagnetism and Heisenberg antiferromagnetism. Comparisons with Monte 

Carlo simulations, RPA result and mean field solutions are also made. 
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I. INTRODUTION 



Low temperature properties of strongly correlated electron systems have been studied 
extensively over the past few decades. The Hubbard model is one of the simplest models of 
interacting electrons, which is believed to capture the main features of electron correlations 
in condensed matter systems. The model is best understood at half hlling, where it exhibits 
the Mott-Hubbard phenomena, namely the existence of insulating behavior accompanied 
by magnetic order in a half-hlled band of electrons. Recent discovery of high-Tc supercon- 
ductivity in doped antiferromagnets (copper oxides) has led Anderson [|^ to suggest that 
superconductivity in these materials is closely related to the physics of strong correlations 
contained in the Hubbard model. The spin-half Heisenberg model (a large U limit of the 
Hubbard model) has formed the basis for a quantitative understanding of the magnetic 
excitations in the insulating stoichiometric phase of the high-Tc materials. Studies of the 
Hubbard and t — J models have also provided qualitative insights into the magnetic and 
transport properties of doped high-Tc materials 0. 

The Hubbard model is dehned by the lattice Hamiltonian 



H = -t Y. Cja + c+ Ci„) +UY - h (na ~h - > (1) 

i i 

where, as usual, cf^ and Cia are the creation and annihilation operators for electrons with 
a 2 ;— component of spin a at lattice site i, and = cf^ Ci„- U is the on site repulsive 
interaction, p the chemical potential, and t the nearest-neighbor hopping amplitude. The 
sum < i,j > is over all pairs of nearest neighbor lattice sites. 

In spite of its apparent simplicity, the Hubbard model is exactly soluble only in one- 
dimension [^. In two-dimensions, relevant to the high-Tc materials, this model has been 
studied by a variety of approximate methods. These include hnite size diagonalization [Q, 
mean- held theory 0, Green’s-function decoupling schemes 0, variational approaches P, the 
random phase approximation (RPA) [^, high temperature series expansion §] and quantum 
Monte Carlo simulations |]9|-[II|. Despite these efforts some of the very basic features of this 
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model, such as its phase diagram in the space of temperature, density and U /t are not fully 
known. The numerical studies do provide a consistent picture for the magnetically ordered 
zero temperature phase at half-hlling. However, away from half-hlling the possibility of long 
range incommensurate magnetic order and superconductivity remain open. 

In this paper we develop a series expansion method for studying the Hubbard model, 
and calculate the ground state energy Eg, local moment L, sublattice magnetization M+, 
the uniform magnetic suscebtibility x± spin stiffness ps, for the half-hlled case, as 
a function of U /t. This study is complimentary to the existing numerical methods and 
provides additional insight into this model. Our results are in broad agreement with the 
Monte Carlo simulations. There are however, some quantitative differences. For example the 
manner in which the local moment L goes to 1/2 as 0, is different in our calculations 
than in the Monte-Carlo simulations. Our results are much closer to the spin-density-wave 
mean field results in this limit. Furthermore, the behavior of magentic susceptibility as a 
function oi U/t changes character at ~ 4 indicating a crossover between spin density 
wave (SDW) antiferromagnet and a Heisenberg antiferromagnet at this value 

The organization of the paper is as follows: In Sec. II we discuss our series expansion 
method for the two-dimensional Hubbard model. In Sec. HI we present results for the 
ground state properties and in Sec. IV for the uniform magnetic susceptibility of the model. 
In Sec. V we derive and compute the spin stiffness of the Hubbard model and estimate the 
spin wave velocity from the relation = Ps/x±- Finally, in Sec. VI, we summarize our 
results. 



II. ISING EXPANSION FOR THE HUBBARD MODEL 

In order to develop series expansion for the Hubbard model, we need to introduce an 
Ising anisotropy into the Hubbard Hamiltonian: 

H, = -Xt Y, (c£c,, + c+ + +/(A) , (2) 

<i,j>,(7 i <i,j> 
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where af = — nj|) is the component of the spin at site i. The particle- hole symmetry 

ensures half-hlling. The function /(A) needs to be chosen such that /(O < A < 1) > 0, 
and /(I) = 0. For A = 0 the atomic limit of the Hubbard model is highly degenerate, 
however the Ising term selects from these the Neel states as the two degenerate ground 
states. Futhermore, this term also introduces a gap in the spectrum at A = 0. For A = 1 
the Ising anisotropy goes to zero and the conventional Hubbard Hamiltonian is recovered. 
Ground state properties of the model for A 7^ 1 can be obtained by an expansion in powers 
of A. If the gap does not close before A = 1 as expected for this model, we can obtain 
properties of the Hubbard model by extrapolating the expansions to A = 1. 

We choose the function /(A) to have the form 

/(A) = J(1 - A) , (3) 

where J is a constant which can be tuned to improve the convergence of the extrapolations. 
In the strong coupling limit, the half-hlled Hubbard model is equivalent to an antiferromag- 
netic Heisenberg model with Hamiltonian: 

H Heisenberg J }j E - 1) , (4) 

<i,j> 

where < i,j > denotes nearest neighbors and Jh = t^ /U. If we choose J in Eq.(3) to be 
much large than Jjj we hnd that the expansion coefficients become very large and oscillate 
rapidly, making extrapolations very difficult. On the other hand, if J is much smaller than 
Jh, the expansion coefficients become very small and a large number of terms will be needed 
to extrapolate to A = 1. The optimum value of J is around Jh- 

In order to employ series expansion techniques, we write the Hamiltonian in Eq.(2) as 

Hy^ = H^ + H' , (5) 

where unperturbed part of Hamiltonian Hq is dehned as 

//o = C/j:(n„-i)(n.j-i) + J ^ tr'y , ( 6 ) 

i ^ ^ <i,j> 
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and the perturbation H' is 



H' = -\t Y. {4aCja + 4^Ci^) - XJ Y 44 ■ (7) 

<i,j>,o <i,j> 

Throughout this paper, we take t = 1 as the energy unit. U and J are two interaction 
parameters. The parameter A simply allows us to extrapolate from the Ising model to the 
half-hlled Hubbard model. 



III. GROUND STATE PROPERTIES 



The ground state energy Eg is dehned as 



B, = < V,(A) I H(\) I V>,(A) > / < V-»(A) I M\) > 

= j;: e„(£/.J)A” , (8) 

n 

where | 'ipg(X) > is the ground state. 

The sublattice magnetization M~^ is dehned as 

2M+ = ^ <Y 4) - ^n) > 

i 

= J) A” , (9) 

n 

where N is the number of lattice sites, e{i) set to be +1 on one sublattice and —1 on the 
other, and the angular brackets refer to ground state expectation values. 

The squared local moment L is dehned via the relation 



L = ^ <Y 

i 

= YUu,J)x^ • 

n 



( 10 ) 



The expansion coefficients for Eg, L, and 2M+ to 11th order in A are listed in tables I-III. 

We hnd that properties which are not sensitive to long range order, such as ground state 
energy Eg and local moment L can be easily obtained to high accuracy with the simplest 
extrapolations. Fade approximants to the function under consideration may be formed 
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directly. In order to estimate the values of Eg and L at particular [/, we compute all Fade 
approximants from series up to 9th, 10th and 11th orders for three different Ising parameter 
J/{t^/U) = 0.75,1.0,1.25, and remove ill-behaved Fade approximants. The estimates are 
set to be the median of the Fade approximants at A = 1, and the errorbars are estimated 
as half the spread of these Fade approximants. 

The values of the ground state energies at A = 1 are Eg = —1.51(2), —1.86(1), —2.934(4) 
and —5.2293(8) for 17 = 1, 4, 10 and 20, respectively. They are plotted in Fig. 1 as solid circles 
with errorbars. The solid line is a guide to the eye. Our calculated ground state energies Eg 
agree with the Monte Carlo simulations. For example, at 17 = 4, the extrapolated ground 
state energy from hnite temperature Monte Carlo simulations is —1.864(1). Frojected 



Monte Carlo simulations |T^ (an algorithm which allows the calculation of ground state 
properties in the canonical ensemble, he., with a hxed number of electrons) gives —1.864(2). 
At 17 = 1, our result {Eg + U/A = —1.26) is close to the slave-boson mean held result 
(^ -1.3) 0 

The squared local moment obtained from the Fade approximants at A = 1 are L = 
0.54(1), 0.75(1), 0.924(2) and 0.9788(4) for 17 = 1, 4, 10 and 20, respectively. They are plotted 
in Fig. 2 as solid circles with errorbars. Once again, the solid line is a guide to the eye. The 
empty squares and dashed lines are the Monte Carlo simulations obtained by White et 
al. 0 at low temperatures [Fig. 6(a) in Ref. 9]. We see that both calculations are in 
good agreement with each other at large U but they show different behaviors at small U. 
Monte Carlo simulations show that local moment L increases rapidly as U increases {concave 
function) while ours indicate a relatively slow growth {convex function ) at small U region, 
as shown in Fig. 3. We would naively expect the series extrapolation to be better at large 
U, and not so well at small U. However, as we now show, our results are rather close to the 
weak coupling results at small 77, which suggest that our convergence is reasonable for all 
U. 

Within the Hartree-Fock approximation, we can derive a simple relationship between 
local moment and sublattice magnetization. The local moment is 
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L =< > + < 7ii > — 2 < >< ni > , 



( 11 ) 



where 



< m > = 0.5 + M+ 

< > = 0.5 - M+ . (12) 

Substituting Eq.(12) into Eq.(ll) we have 

L = 2(M+)2 + 0.5 , (13) 

The mean held solution for the sublattice magnetization is |^,|T^ 

M+ = A/U , (14) 



where A can be obtained by solving the gap equation 



U 



E 



1 



2N ^ {el + A^y/^ 



= 1 



(15) 



N is the total number of sites and e(k) = —2t{coskx + cosky). The local moment given 
by numerical solution of the gap equation is plotted in Fig. 2, a dashed line. One can see 
that the function L{U) has a convex shape at small U which agrees with our conclusion. At 
U = 1, our result {L = 0.54) agrees with the slave-boson mean held result {L ^ 0.53) H 
The extrapolation of sublattice magnetization M+ appears much harder. In the 2D 
spin-1 square lattice Heisenberg antiferromagnet, the sublattice magnetization (2M+) is 
estimated by expansions around the Ising limit. In order to extrapolate the 2M+ series 
reliably, one has to remove the singularity of the form [1 — {J±/J\\)'^]^ caused by Goldstone 
modes HI. where J± is the exchange perpendicular to the direction of ordering and J\\ the 
exchange along it. One way of doing this is to go a new variable 6 given by 0 






( 16 ) 



We expect to have the same situation for the sublattice magnetization 2M^ in the 2D 
Hubbard model because it is equivalent to the Heisenberg model at large U. Here we set a 
new varible 6 to be 
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( 17 ) 



where 5 = 1 at A = 1. Then the series for the sublattice magnetization 2M+ becomes 
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2M* = Y, bn{U,J)y' 



( 18 ) 



n=0 



where bn{U, J) can be calculated from previous series mn{U, J). Again, we compute all Fade 
approximates from 7th order to 11th order series. The median of Fade approximants at 
5 = 1 give 2M+ = 0.12(5), 0.37(5), 0.58(3), and 0.61(3) for 7/ = 1,4, 10 and 20, respectively. 
They are plotted in Fig. 3, solid circles with errorbars. The solid line is a guide to the eye. 

At large 7/, the effective Hamiltonian for the half-hlled Hubbard model is known to be 
the spin-A Heisenberg model. The quantum zero-point spin wave fluctuations reduce the 
sublattice magnetization of 2M+ from the classical value of unity to about 2M+ = 0.6 as 



hrst obtained by Anderson [18]. This is also supported by Monte Carlo simulations [19] and 



series expansion studies [^. At 77 = 20, we estimate 2M+ = 0.61 which is close to the value 
0.6 of the Heisenberg model. Also, at 7/ = 4, we can compare our results [2M+ = 0.37(5)] 
with the Monto Carlo simulations [2M+ = 0.39(5)] of Hirsch and Tang and hnd fairly 
good agreement. 

Also, we compare our results with the mean held solutions and the RFA calculations. 
The sublattice magnetization given by numerical solution of the gap equation is plotted in 
Fig. 3, as a long dashed line. At large 77, the mean held solution gives 2M+ = 1, which is 
far from the true value (2M+ = 0.6) of the Heisenberg model because of the neglect of the 
quantum huctuations. In order to take into account the huctuation ehects on the sublattice 
magnetization, one simply has to calculate the self-energy correction to the one particle 
Green’s function as done by Schrieher, Wen and Zhang in random phase approximation. 
Their RFA result is plotted in Fig. 3 by a dotted line. Clearly, the RFA result brings the 
mean held result much closer to our numerical estimates. 



We note again that the purpose of introducing the Ising anisotropy J(1 — A)crf ct| into 
the conventional Hubbard Hamiltonian is to be able to carry out the series expansions. The 



ground state properties of the Hubbard model should not depend on the parameter J of the 
Ising anisotropy. One can also choose another type of function /(A) in Eq.(2): 

/(A) = J(1-A^) . (19) 

This form of /(A) leads to an expansion in powers of A^. This allows one to treat the J 
term generated by the pertnrbation in the Hubbard model and the Ising term introdnced 
by hand on an equal footing. However, it has the disadvantage that the number of series 
coefficients available for extrapolation is rednced by half. The resnlting extrapolations at 
A = 1 are consistent with those obtained earlier. 

IV. MAGNETIC SUSCEPTIBILITY 

We now discuss the behavior of the uniform magnetic susceptibility of the 2D Hub- 
bard model from the series expansions. We tnrn on the magnetic held and add a Zeeman 
interaction into the 2D Hubbard Hamiltonian: 

Hxh = -\t Y. (4a Cja + C+, Cia) + U^ 

<ij>,a i ^ ^ 

+ J(1-A) i: - hJ2af , (20) 

<i,j> i 

where h is magnetic held and is the x— component of the spin at site i. The nniform 
magnetic susceptibility yj_ is dehned as 

4y± = -d^E/dh^ \h=o 

= Y^n{U,J)\^ . ( 21 ) 

n 

The expansion coefficients of A\± to 9th order in A are given in Table IV. 

At half-hlling, 2D Hnbbard model has antiferromagnetic long range order. One can see 
that x± series is dominated by a simple pole at A = — 1. We need to remove this singnlarity 
before fnrther analysis. This is done by going to a new variable f3 given by 

/3 = 2A/(1 + A) . (22) 
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The resulting series has the form 



4x± = E 2/n(f/,J)/S" , (23) 

n=0 

where yn{U, J) can be obtained from the previous series J). But this new series still 

contains singularities at the Hubbard point (A = 1 or /3 = 1). Going to a new variable 5 
dehned by 

1_5=(1_/3)1/2 ^ (24) 



and constructing Bade approximants, we hnd that only two Fade approximants {K = 4, M = 
5 and K = 5, M = 4) from 9th order series show good behavior and both give the same 
results. For this reason it is difficult to estimate errors in these resnlts. The nniform magnetic 
susceptibility obtained from these Fade approximants are 4x± = 0.94, 0.58, 0.75, 0.92 and 
1.40 for U = 1,4,6, 10 and 20, respectively. They are plotted in Fig. 4 as solid circles. The 
solid line is a guide to the eye. 

At U = 20, x±/U = X±J = 0.0175 is close to the large U limit value x± Jh = 0.0163 
of the Heisenberg model shown by a dashed line in Fig. 4(a). At H = 4, onr resnlt 
{4x1. = 0.58) agrees with the estimated valne (~ 0.53) of Monte Carlo simulations at low 



temperatnre by White et al. |^, bnt does not agree with the estimated value (~ 1.0) of 



projected Monte Carlo simulations by Furukawa and Imada |^. At U = 10, onr resnlt 
{4xi. = 0.92) agrees with the estimated valne (~ 0.98) of Monte Carlo simulations by 
Moreo |2^ . 



One can see that onr magnetic susceptibility data indicate a dip at f/ = f/c ~ 4, as shown 
in the Fig. 4(b). It suggests that there is a crossover around Um in the behavior of the 2D 
Hubbard model at half filling. For U > Uc^ we can see that magnetic susceptibility {4xi_) is 
almost proportional to U or inversely proportional to the spin snperexchange J . This can be 
taken as evidence that the magnetic state for U > Uc for 2D Hubbard model is that of the 
Heisenberg model. Bnt, for U < Uc, magnetic susceptibility {4xi_) has a different behavior, 
and decreases as U increases. This can be interpreted as a weak conpling behavior of the 
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SDW ground state. Thus we obtain a crossover between these types of behavior around 
U ^ A. 

Also, we compare our result with the mean held solutions of uniform magnetic suscepti- 
bility [0: 

Xi = T[(1-AW/,-(7/|//i)-‘;-1| . (26) 

where 

I = — V ^ 

^ 2N ^ E^{k) 

r ^ L V 

^ 2N Y ^^(k) ■ 

Here T^(k) = [e^(k) -|- It is plotted in Fig. 4(b), a solid line. We see that the mean 

held result is qualitatively correct, but overestimates x±^ especially at large U. 



V. SPIN STIFFNESS OF THE HUBBARD MODEL 



The spin stihness constant ps is a measure of the response of the spin system in an 
ordered phase to a twist of the order parameter. To our knowledge, there are no quan- 
tum Monte Carlo estimates for the spin stihness in the Hubbard model. Derivation and 
calculation of this quantity for the Hubbard model was discussed by us recently in a short 
communication For completeness, we discuss it again here. 

If we rotate the ordering direction by a small angle 6 along a given direction such as y 
axis, then the spin stihness constant ps can be dehned through the increase of the ground 
state energy: 



E,(») = E,(» = 0) + + O{0^) , 



(26) 



This rotation can be carried out by the following transformation applied to the Fermion 
operators: 

/ \ / 

cos(f)/2 sin(j)/2 
—sin 0/2 cos(!)/2 



( 


( 








V 



\ 




/ 





( 27 ) 
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and 



fcr] 


/ 






Vi J 


V 



cos ^12 sin(f)/2 
—sin 0/2 cos(f)/2 

This transformation generates desired spin rotation: 



\ 


fd) 


/ 


di" ! 



c'^c[ 

= Gz COS (f) + Ux sin 0 
ax' = c/" c'l — c'/“ c'l 

= —az sin (f) + ax cos 0 



( 28 ) 



( 29 ) 



After rotation by a relative angle 6 ( that is letting 0 change hy 9/2 ) between neighboring 
sites separated along y axis {y denotes unit distance in y direction) , H in Eq.(l) becomes 

He = H + + 0{9^) , ( 30 ) 



where 



x02 

Ljdia — („+ \ 

o '^t+ycT W f^i+ya- ) 

i,a 

i 



( 31 ) 



The “diamagnetic” term is already of order 6^ so for the calculation of the energy 

to order 0^, it can be replaced by its expectation value in the ground state of the 9 = 0 
Hamiltonian. We have 



P 



dia 

s 



1 

8 

1 

8 



{Kinetic Energy) 

[E,(» = 0) - I (n - L)| 



(32) 



where n is a band filling, and L a local moment defined in Eq.(lO). We have already obtained 
Eg and L in Sec. Ill, and n = 1 for half-filling. 
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The contribution of the “paramagnetic” term to the ground state energy in order 



can be obtained by treating it in second order perturbation theory. We have 

d'^E 

nara _ oy_“ I 

^ QQ2 1^-0 ’ 

= Y.PnX' , 



(33) 



where E represents energy of the Hamiltonian E[\ + The series coefficients of 

to 9th order in A are given in the Table V. 

We use the same Fade analysis for as for uniform magnetic susceptibility. We obtain 
the spin stiffness ps = pf°‘ + = 0.186(15), 0.15(1), 0.077(7) and 0.039(5) for 7/ = 1, 4, 10 

and 20, respectively. They are plotted in Fig. 5 as hlled circles with a dashed line as guide 
to the eye. The errorbars represent the spread in the Fade estimates. At f/ = 20, our result 



{psU = 0.78) agrees well with the value (0.73) of the Heisenberg model shown by a short 
solid line. 

We also compare our results with the Hartree-Fock approximation for the spin stiff- 
ness 



21 



Ps = 



t ^ e(k) cos 



IN y £(k) 



(34) 



The spin stiffness given by mean field solution of the gap equation is plotted in Fig. 5(b) as 
a solid line. One can see that the mean field approximation overestimates the stiffness at 
large U . AiU = 20, the mean field result is psU = 0.98 compared to 0.73 for the Heisenberg 
model. 

We note that the spin stiffness (filled squares in Fig. 1) from the variational Monte 
Carlo method with a Gutzwiller-type wave function are even larger than values of the mean 



field solution |21]. We believe that the large discrepancy is due to the missing spin flip 



processes in the Gutzwiller variational wave function used in |^. Their calculations only 
get contributions to the spin-stiffness beyond the Hartree-Fock result from the “diamagnetic” 
term, whereas the “paramagnetic” part of the spin stiffness which contains spin-flip processes 
does not get corrected. The expectation value of the operator between the Hartree- 



Fock and Gutzwiller wave functions is zero [^. This suggests that one might improve the 
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variational Monte Carlo method by adding a spin flip operator Ogp into the Gutzwiller-type 
wave function: 



'^GW > = H [1 “ (1 “ + P ^ Osp] I > 

i 



where p is a parameter and 



Osp (Cj| Ci^yi 

We now turn to the calculation of spin-wave velocity. This can be calculated analytically 
by adding RPA fluctuations to the SDW state [0. Designed, at hrst sight, for intermediate 
and weak coupling {U/At < 1) this method also interpolates smoothly to the Heisenberg 
limit at large U/t. 

Alternatively, there is a relationship between the spin stiffness Ps of the order param- 
eter with respect to a spiral twist, the spin wave velocity Vs and the uniform magnetic 
susceptibility xj_ transverse to the stagged spin density: 

Ps =vjx± ■ (35) 



It was originally derived by the hydrodynamic theory of Halperin and Hohenberg . The 
spin wave velocity obtained from this relationship (p^ = 'VgX±) ci-re Vg = 0.89, 1.02, 0.58 and 
0.34 for 17 = 1, 4, 10 and 20, respectively. They are plotted in Fig. 6 as hlled circles with 
a solid line as a guide to the eye. At 17 = 20, our Vg = 0.34 agrees well with the limiting 
value {vg = 1.18^/2J' = 0.33) of the Heisenberg model |l^, a dashed line in Fig. 6, where 
J' = AtyU. 

We compare this result with the RPA solution of Schrieffer, Wen and Zhang [0 : 



V, = - A^U h) , 



where 



^ 1 (kx) 



2N ^ £3(k) 



( 36 ) 



( 37 ) 



It is plotted in Fig. 6 as a dotted line. We see that the RPA result for the spin wave velocity 
are lower than our series expansion results. 
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VI. SUMMARY 



To summarize, in this paper we developed series expansions for the half-hlled Hub- 
bard model, by introducing an Ising anisotropy into the conventional Hubbard Hamiltonian. 
Series were developed for the ground state energy, the local moment, the sublattice magne- 
tization, the uniform magnetic susceptibility and the spin stiffness, and were extrapolted to 
the Hubbard model using standard series extrapolation methods. In general, our results are 
in good agreement with the quantum Monte Carlo simulations. There are some quantitative 
differences, which are also discussed. Our magentic susceptibility data indicate a crossover 
between SDW antiferromagnetic state at small U and a Heisenberg antiferromagnet at large 
U at U ~ 4. We note that our results for ps as well as Ml are somewhat higher than the 
known results for the Heisenberg models. We believe this reflects the fact that the series 
have not converged as well as for the Heisenberg model and hence the reduction in these 
quantities due to the zero-point spin-wave fluctuations is not fully accounted for. Still, the 
convergence (~ 5%) is quite reasonable. 
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TABLES 



TABLE I. Expansion coefficients for the ground state energy. 



u 


JKfju) 


en(U, J), n = 0 — 11 


1 


0.75 


-3.25 3.0 -0.64 -0.5376 -0.36882596 -0.17599739 -0.0037101999 0.11731865 
0.17602728 0.17413953 0.11347540 -0.0088559404 


1 


1.0 


-4.25 4.0 -0.5 -0.4375 -0.34273695 -0.23157505 -0.12102132 -0.025358881 
0.046126082 0.089325558 0.10390817 0.091603525 


1 


1.25 


-5.25 5.0 -0.41025641 -0.36817883 -0.30808094 -0.23717612 -0.16320812 
-0.093237945 -0.032792472 0.014501216 0.046805520 0.063642701 


4 


0.75 


-1.75 0.75 -0.75294118 -0.18602076 0.032575733 -0.11521776 0.051494035 
0.41270231 0.36649252 -0.23825557 -0.77368699 -0.75051947 


4 


1.0 


-2.0 1.0 -0.69565217 -0.21172023 -0.030847663 -0.060114049 -0.011016385 
0.16603971 0.22775615 0.035681298 -0.22467766 -0.32023724 


4 


1.25 


-2.25 1.25 -0.64646464 -0.22854811 0.042172811 -0.035358477 -0.023047559 
0.066640645 0.13361749 0.080005270 -0.048764273 -0.13678257 


10 


0.75 


-2.8 0.3 -0.38004751 -0.018957239 0.066441099 -0.087012637 0.038082237 
0.16089588 0.11974079 -0.075315493 -0.24874869 -0.20497804 


10 


1.0 


-2.9 0.4 -0.37383178 -0.024456284 -0.045510051 -0.064243043 0.0053501507 
0.075910665 0.074634170 0.0084898042 -0.068321715 -0.096846836 


10 


1.25 


-3.0 0.5 -0.36781609 -0.029594398 0.033584424 -0.050422135 -0.0070731834 
0.039086550 0.046846607 0.021090155 -0.015763853 -0.040974167 


20 


0.75 


-5.15 0.15 -0.19740901 -0.0025574270 0.041371562 -0.044384097 0.016313799 
0.075735746 0.056531614 -0.041034781 -0.11826259 -0.071871255 


20 


1.0 


-5.2 0.2 -0.19656020 -0.0033806422 -0.030305728 -0.033253249 0.0010839733 
0.034576709 0.03478823 0.019973437 -0.033086550 -0.037873386 


20 


1.25 


-5.25 0.25 -0.19571865 -0.0041896960 -0.023690877 -0.02656723 -0.0044976664 
0.016991751 0.021350382 0.0088434956 -0.0081635325 -0.016372525 
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TABLE II. Expansion coefficients for the squared local moment. 



u 


JKflv) 


ln{U, J), n = 0 — 11 


1 


0.75 


1.0 0.0 - 0.2048 - 0.34406400 0 .34363876 - 0.18584371 0.073558744 
0.3 2646019 0.45962201 0.39942169 0.126 77212 - 0.33704683 


1 


1.0 


1.0 0.0 - 0.125 - 0.21875 - 0.25342785 - 0.21761601 - 0.12058227 
0.010355528 0.13808712 0.22622422 0.248 12902 0.19219824 


1 


1.25 


1.0 0.0 - 0.084155161 - 0.15104772 0 .18802143 - 0.18850833 - 0.15340922 
- 0.090723264 - 0.013648193 0.062253229 0.12199406 0.15344711 


4 


0.75 


1.0 0.0 - 0.28346021 - 0.14006269 0.020099625 - 0.042208513 0.13810651 
0.6 6290690 0.47876441 - 0.91546939 - 1.99 197718 - 1.29230064 


4 


1.0 


1.0 0.0 - 0.24196597 - 0.147283644 0.0033966659 0.0020465901 0.041625973 
0.27003095 0.33545751 - 0.10812405 - 0.68716227 - 0.74045463 


4 


1.25 


1.0 0.0 - 0.20895827 - 0.14774827 0.017014379 0.014267802 0.020573983 
0.1 1999820 0.20351777 0.071247345 - 0.22059729 - 0.38236644 


10 


0.75 


1.0 0.0 - 0.072218053 - 0.0072046514 0.022443265 - 0.034066355 0.041464469 
0.1 1664774 0.058849733 - 0.12141417 - 0.25211532 - 0.14201779 


10 


1.0 


1.0 0.0 - 0.069875098 - 0.0091425362 - 0.014415554 - 0.024755587 0.016845440 
0.060456748 0.045459350 - 0.021365905 - 0.086403185 - 0.092929358 


10 


1.25 


1.0 0.0 - 0.067644339 - 0.010885296 0.0099895925 - 0.019054664 0.0064590300 
0.035218951 0.032572341 0.0029017356 - 0.031199252 - 0.048202015 


20 


0.75 


1.0 0.0 - 0.019485158 - 0.00050485912 0.0079729731 - 0.0088647321 0.0094722236 
0.027330357 0.013712645 - 0.031211772 - 0.058546614 - 0.016716436 


20 


1.0 


1.0 0.0 - 0.019317955 - 0.00066449970 - 0.0057661935 - 0.0066346164 0.0037510047 
0.013829393 0.010568695 - 0.0058676230 - 0.020309374 - 0.016500150 


20 


1.25 


1.0 0.0 - 0.019152896 - 0.00082000166 - 0.0044501566 - 0.0052935510 0.013853661 
0.078617712 0.0074609205 0.00033154323 - 0.0074336626 - 0.0092115763 
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TABLE III. Expansion coefficients for the sublattice magnetization. 



u 


JKflv) 


mn{U, J), n = 0 — 11 


1 


0.75 


1.0 0.0 -0.2048 -0.34406400 0.37832834 -0.30858617 -0.16797349 0.0039694233 
0.17811940 0.32196832 0.35932350 0.14690771 


1 


1.0 


1.0 0.0 -0.125 -0.21875 -0.26577581 -0.26255537 -0.21639940 -0.14085987 
-0.050839444 0.039513722 0.11520618 0.15738379 


1 


1.25 


1.0 0.0 -0.084155161 -0.15104772 .19346230 -0.20866262 -0.19818728 -0.16672332 
-0.12065102 -0.066877941 -0.012424680 0.035149309 


4 


0.75 


1.0 0.0 -0.28346021 -0.14006269 0.29094791 -0.70653657 -0.11899194 1.78106763 
2.65518980 0.17150140 -4.57860227 -8.08441904 


4 


1.0 


1.0 0.0 -0.24196597 -0.14728364 -0.16340074 -0.36296753 -0.23038035 
0.49230029 1.09523777 0.69655017 -0.61481862 -2.00113858 


4 


1.25 


1.0 0.0 -0.20895827 -0.14774827 0.11819956 -0.21655277 -0.20123752 0.10926874 
0.47800300 0.49583051 0.078098871 -0.50046276 


10 


0.75 


1.0 0.0 -0.072218053 -0.0072046514 0.38636976 -0.78484180 -0.10195808 
1.63572360 2.29944530 -0.28902064 -4.69346784 -5.25870098 


10 


1.0 


1.0 0.0 -0.069875098 -0.0091425362 -0.21519753 -0.44065350 -0.20500969 
0.47694043 0.93375402 0.48462349 -0.74441315 -1.65106094 


10 


1.25 


1.0 0.0 -0.067644339 -0.010885296 0.13638406 -0.28150839 -0.18799630 
0.13527913 0.41216832 0.36759744 -0.024358458 -0.48058650 


20 


0.75 


1.0 0.0 -0.019485158 -0.00050485912 0.39342412 -0.78866044 -0.12416041 
1.59373845 2.17733968 -0.58972155 -4.77779001 -4.11947852 


20 


1.0 


1.0 0.0 -0.019317955 -0.00066449970 -0.22095641 -0.44335646 -0.21769894 
0.45233967 0.87748571 0.35719439 -0.84087768 -1.42111547 


20 


1.25 


1.0 0.0 -0.019152860 -0.00082000166 -0.14118233 -0.28358286 -0.19603500 
0.11902603 0.38038377 0.30147500 -0.092194212 -0.43888627 




TABLE IV. Expansion coefficients for the uniform magnetic susceptibility. 


U 


JliflU) 


Xn(U, J), n = 0 — 9 


1 


1.0 


0.25 0.25 0.17261905 0.027848648 -0.1262274 -0.19962669 -0.11626869 
0.12849035 0.43255914 0.60788905 


4 


1.0 


1.0 1.0 -0.67116442 -2.82447457 -1.58009047 4.43513119 8.64219154 
1.42302418 -15.39992391 -22.54729517 


6 


1.0 


1.5 1.5 -1.48149027 -4.91068070 -1.92590553 8.56642180 13.73914475 
-1.78852254 -29.87404955 -32.12999668 


10 


1.0 


2.5 2.5 -2.99306647 -8.80889678 -2.67462512 16.12944762 23.047960372 
-8.54214685 -56.78749173 -44.88013932 


20 


1.0 


5.0 5.0 -6.48980810 -18.15418203 -4.79124584 33.98201522 45.62214165 
-23.66505149 -119.66330645 -76.12120519 
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TABLE V. Expansion coefficients for the “paramagnetic” part of the spin stiffness. 



u 


JliflU) 


Pn{U, J), n = 0 — 9 


1 


1.0 


0.0 0.0 0.0 0.0 0.000047115128 0.00023699734 0.000655001 
0.0012665331 0.0018294204 0.0019236148 


4 


1.0 


0.0 0.0 0.0 0.0 0.0089710818 0.028888083 0.029891453 
- 0.020672065 - 0.090109640 - 0.087959801 


10 


1.0 


0.0 0.0 0.0 0.0 0.012111686 0.027777667 0.018651461 
- 0.018905322 - 0.049016789 - 0.033652011 


20 


1.0 


0.0 0.0 0.0 0.0 0.0076518484 0.015940557 0.010548308 
- 0.0087795777 - 0.023213629 - 0.013514588 



21 





Figure Captions 



Fig. 1 The ground state energy as a function of U. Filled circles are our calcu- 
lated data, with a solid line as guide to the eye. Errorbars indicate spread of the Fade 

approximants. Here Eg + U /A = Kinetic Energy + U < n^ ni >. 

Fig. 2 The squared local moment L as a function of U. Filled circles are our cal- 
culated data, with a solid line as guide to the eye. Errorbars indicate spread of the Fade 



approximants. Open squares are the results of Monte Carlo simulation by White et al. [10 



with dotted line as guide to the eye. The dashed line is the result of the mean held solution. 

Fig. 3 Sub lattice magnetization 2M+ as a function of U. Filled circles are our 
calculated data, with a solid line as guide to the eye. Errorbars indicate spread of the Fade 
approximants. Long dashed line and dotted line represent mean held solution and RFA 
result respectively. A short dashed line indicates the limiting large U (Heisenberg) value. 

Fig. 4 (a) The ratio of the uniform magnetic susceptibility to f/ as a function of 

U. Filled circles are our calculated data, with a dashed line as guide to the eye. A short 
dashed line indicates the limiting large U (Heisenberg) value, (b) The uniform magnetic 
susceptibility as a function of U. Filled circles are our calculated data, with a dashed line 
as guide to the eye. A crossover between SDW antiferromagnetic state and Heisenberg 
antiferromagnet is around f/c ~ 4. 

Fig. 5 (a) The product of the spin stihness and f/ as a function of U. Filled circles 

are our calculated data, with a solid line as guide to the eye. A short solid line indicates the 
limiting large U (Heisenberg) value, (b) The spin stihness as a function of U. Filled circles 
are our calculated data, with a solid line as guide to the eye. Errorbars indicate spread of 
the Fade approximants. The solid line is the mean held result. The hlled squares are results 
of variational Monte Carlo simulations . 

Fig. 6 The spin wave velocity as a function of U. Filled circles are our calculated 
data, with a solid line as a guide to the eye. The dotted line and dashed line represent 
results of RFA and Heisenberg model respectively. 
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